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ABSTRACT 

The kinematics of damped Lyman-a absorbers (DLAs) are difficult to reproduce in hi- 
erarchical galaxy formation models, particularly the preponderance of wide systems. 
We investigate DLA kinematics at z — 3 using high-resolution cosmological hydro- 
dynamical simulations that include a heuristic model for galactic outflows. Without 
outflows, our simulations fail to yield enough wide DLAs, as in previous studies. With 
outflows, predicted DLA kinematics are in much better agreement with observations. 
Comparing two outflow models, we find that a model based on momentum-driven 
wind scalings provides the best match to the observed DLA kinematic statistics of 
Prochaska & Wolfe. In this model, DLAs typically arise a few kpc away from galaxies 
that would be identified in emission. Narrow DLAs can arise from any halo and galaxy 
mass, but wide ones only arise in halos with mass > 10 11 Mq, from either large central 
or small satellite galaxies. This implies that the success of this outflow model origi- 
nates from being most efficient at pushing gas out from small satellite galaxies living 
in larger halos. This increases the cross-section for large halos relative to smaller ones, 
thereby yielding wider kinematics. Our simulations do not include radiative transfer 
effects or detailed metal tracking, and outflows are modeled heuristically, but they 
strongly suggest that galactic outflows are central to understanding DLA kinemat- 
ics. An interesting consequence is that DLA kinematics may place constraints on the 
nature and efficiency of gas ejection from high-z galaxies. 
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1 1 INTRODUCTION 

Damped Lyman alpha systems (DLAs), i.e. Hi absorp- 
tion line syste ms having colum n densities of Nm > 
2 x 10 20 cm" 2 (|Wolfe et al. I Il986f ), contain the dominant 
reservoir of neutral hydrogen in the Univ erse (see re- 
view by IWolfe. Gawiser. fc Prochaska I l2005h . Since neu- 
tral hydrogen is closely connected with star formation (e.g. 

iKennicutt | |l998h . DLAs represent the neutral gas repos- 
itory for fuelling star formati on. Influential work by 

IStorrie-Lombardi fe Wolfel (2000) suggested that the global 
neutral gas content measured from DLAs declines with 
cosmic epoch in concert with the growth in cosmic stel- 
lar mass. Other studies with different selection meth- 
ods s uggest that the decline in neu tral gas is not as 
steep jRao. Turnshek. fc Nestor"! [2006'). Nevertheless, one 
expects some connection between galaxies identified in stel- 
lar emission and those identified in gas absorption such as 



DLAs. Studying this connection has been the focus of a large 
number of investigations, enabled by the ability to routinely 
catalogue absorbers and galaxies at z > 2 where optical 
Lya absorption line studies have accurately characterised 
the DLA population. 

Despite these efforts, a deep understanding of the 
relationship between DLAs and emission-selected galax- 
ies remains elusive. Studies of neutral hydrogen ab- 
sorption around Lyman break-selected galaxies (LBGs) 
at z ~ 2 — 3 show a wide range of absorption 
strengths (|Adelberger et al. I l2005h . DLAs have gas-phase 
metallicities that are lower than that seen in LBG s, and 
show a much larger spread (|Prochaska et al~l I2005T ) . Mea- 
sures of the C n* cooling rates in DLAs indicate that the 
star form ation rates are dra matically lower than that seen 
in LBGs (|Wolfe et al. Il2004 ). Low-redshift imaging of DLAs 
show a heterogeneous population of generally sub-i* galax- 



© 0000 RAS 



2 Hong et al. 



ies, but it is unclear if this is relevant to high- z DLAs. On the 
other hand, the clustering of LBGs an d DLAs at z ~ 2 — 3 
suggests that they occupy similar halos (jBouche et al. Il2005l : 
ICooke et al."1l2006h . Hence, even though a wealth of obser- 
vational data exists for DLAs, it cannot be definitively said 
whether DLAs are lower-mass LBG analogs, whether they 
are a phase of galaxy evolution that eventually leads to 
LBGs, or whether they are an altogether separate galaxy 
population. 

One approach for constructing a unified framework for 
the nature of DLAs is to employ numerical simulations that 
account for the hierarchical growth of galaxies including gas 
physical processes a nd star formation. Early investigations 
using simulations bv lHaehnelt. Steinmetz. k, RauchTl|l998l . 
hereafter HSR98) suggested that DLAs traced clumps of 
cold gas that were in the process of assembling into a larger 
galaxies. This proto-galactic clump (PGC) model favoured 
the interpretation of DLAs as a precursor phase to LBGs. 
iGardner et aTl |200ll ) used simulations in cosmological vol- 
umes to argue that the bulk of DLAs must arise near dwarf 
galaxies, suggesting that they are low-mass LBG analogs. 
iMaller et al7T(|200ll ) used semi-analytic models to argue that 
DLAs must come from gas distributions within dark ha- 
los that are extended relative to that expected from an- 
gular momentum support, and hence another mechanism 
such as tidal stripping may be important for obtaining the 
correc t cross-section for DLA absorption. iNagamine et al.l 
120071 ) found that DLA properties are significantly affected 
by galactic outflows, and argued that the more successful 
models showed DLAs as being lower- m ass analogs of LBGs. 
Simulations bv lRazoumov et al. (2006) incorporated radia- 
tive transfer to model the neutral gas distribution more ac- 
curately, and found that DLAs can arise from a variety of en- 
vironments, from the centres of small galaxies to filam entary 
tidal structure in large halos. iPontzen et al. I ((2008) used 
hydrodynamic simulations with star formation and outflows 
together with radiative transfer, and were able to match 
the column densities and (for the first time) the metallicity 
range of DLAs. While these successes are impressive, one set 
of observations have consistently confounded all such mod- 
els: The kinematic properties of DLAs. 

IProchaska fc Wolfe I (|l997l . hereafter PW97) observed 
the detailed kinematics of DLAs via their low-ionization 
metal lines. They devised four kinematic diagnostics to 
quantify their findings, the key one being the velocity ex- 
tent of DLAs. The observations show a pronounced tail to 
large velocity extents. Such a tail is not seen in simu la- 
tions, from the earliest models llProchaska fc Wolfe ll200ll) to 
the most sophisticate d recent ones i Razoumov et al. 1120081 ; 
IPontzen et al. I I2OO8I ). In general, to match the observed 
kinematics the absorption cross section of low-mas s halo s 
must be reduced significantly (|Barnes fc HaehneltJ [2009). 
which does not arise naturally in fully dynamical models. 

Fundamentally, the difficulty with all these models is 
that they yield few absorbers whose velocity spread signifi- 
cantly exceeds the characteristic velocity of the halo in which 
the system resides. Since early gas-rich galaxies tend to be 
small, this makes it difficult to reproduce the observed high- 
velocity (> 200 km/s) tail in DLA system widths. This im- 
plies that the gravitational dynamics of hierarchical assem- 
bly alone cannot explain the distribution of the observed 
kinematics. Hence, any model that does not add a signifi- 



cant component of non-gravitational velocities to the neutral 
gas seems doomed to fail the DLA kinemat i cs test . Alter- 
natively, the original model of Wolfe et al. (1986), which 
proposed that DLAs are large, puffy rotating disks, pro- 
vides an excellent agreement with the observed kinematics, 
but forming enough such objects by z ~ 3 is challenging in 
cu rrently-favoured cosmologie s . Hen ce the oft-debated claim 
by IProchaska fc Wolfe 1 (|l997l . l200ll ) that cold dark matter 
(CDM) models cannot straightforwardly reproduce the kine- 
matic structure of DLAs remains viable. 

While it could be that the currently-favoured hierar- 
chical galaxy formation scenarios are incorrect, their wide- 
ranging successes in other areas suggest that the explanation 
probably lies somewhere within the poorly-understood bary- 
onic physics associated wi th galaxy formation. For DLAs, it 
has been suggested (e.g. IPontzen et al. I [20081 ) that galac- 
tic outflows may be the missing ingredient needed to in- 
crease the kinematic widths and to lower the absorption 
cross-sections of smaller halos. Recent work indicates that 
essentially all star-forming galaxies at z > 1 are ge n- 



erating outflows of hundreds of km/s (e.g . iPettini et al 



2001; ISteidel et al.ll20oi IShaplev et al.ll2005l ; IWeiner et a" 



20091) . with mass outflow rat es comparable to or more than 



their star formation rates l|Erb et alj 120061 : ISteidel et alj 
2010). Theoretically, outflows are bel ieved to be responsible 
for sup pressing early star formation (|Springel fc Hernquist I 
l2003bf ). without which the cosmic stellar mas s density would 
far e xceed observations (|Dave et al. I l200ll ; iBalogh et al. I 
120011 ). But it is unclear how such outflows affect DLAs. Their 
star formation surface densities are quite low (| Wolfe et al. I 
|2004 ). so they may not generate outflows at all. Even if they 
do, it is unclear how much mass is being carried in outflows, 
and whether it would result in enhanced DLA kinematic 
widths in accord with observations. Simulating the impact 
of outflows on DLA kinematics is a challenging problem, be- 
cause it requires both redistributing neutral gas accurately 
as well as self-consistently modeling the energy deposition 
by outflows into the surrounding gas. 

Recently, several groups have considered the impact 
of outflows on DLA pr operties, specifically kinematics. 
iBarnes fc HaehneltJ (2009) used a semi-analytic model to in- 
terpret observations of DLAs and faint extended Lya emit- 
ters, and showed that the incidence rate and kinematics of 
DLAs are best reproduced if neutral gas is removed from 
halos with circular velocitie s < 50 — 70 km/s, s uggesting an 
important r ole for outflows. | Tescari et al.l (|2009l ) followed up 
the work of INagamine et aL (2007) to assess, among other 
things, the impact of various outflow models on DLA kine- 
matics, and found that none of their outflow models were 
able to reproduce the kinematics completely, although a 
mo del based on momen tum-driven wind scalings came clos- 
est. IZwaan et al.l (|2008l ) argued from an observational per- 
spective that DLA kinematics must be influenced by out- 
flows, by comparing the Hi kinematic widths of local Hi 
galaxies and those of high-z DLAs, and showing that (as- 
suming that locally-identified DLA galaxies are representa- 
tive of those at high-z) some other process such as outflows 
must be enhancing the kinematics of high-z DLAs. Hence 
there is growing evidence that DLA kinematics are strongly 
impacted by galactic outflows. 

In this paper, we investigate whether galactic outflows 
can both lower the absorption cross-section of small halos 
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and broaden the velocity widths enough to bring the the- 
oretical DLA kinematics into accord with the observations. 
To do so, we employ cosmological hydrodynamic simulations 
of galaxy formation, with the key addition being several 
heuristic models of galactic feedback. We consider three vari- 
ants: No outflows; a model with constant outflow speed and 
mass loading factor (i.e. the mass outflow rate rela tive to the 
star formation rate; ISpringel fc Hernquist |[2003bl . hereafter 
SH03); and a model where the outflow speed and mass load- 
ing factor sca le according to expectations from mome ntum- 
driven winds (|Murrav. Quataert. fc Thompson~| [2005). Such 
momentum-driven wind scalings implemented within our 
models have proved remarkably successful at reproducing 
a wide range of observations, includ ing: the observed IGM 
meta l content fr om z ~ — 6 ijOppenheimer fc Dave I 
20061. |200S| . |2009| ); observ ations of early galaxy luminosity 
funct i ons and evolution |Pave. Finlator. fc Oppenheimer I 
I2006L n ouwens et al. 1 I2007T): the g alaxy mass-metallicity 
relationship ( Finlator fc Pavel 120081); the observed enrich- 
ment in various baryonic phases l|Pave fc Oppenheimer I 
t he metal and entropy content of low-r edshift galaxy 
groups |Pave. Oppenheimer. fc Sivanandam II2008T) ; th e H I 
distribution in the local universe ((Popping et al. 1120091 ) : and 
the f aint-end slope of the pres ent day stellar mass func- 
tion (|Oppenheimer et all 120101 ). Momentum-driven wind 
scalings of out flow velocity are also observed in n earby star- 
burst outflows ( Martin 2005; Rupkc et al. 2005), providing 
an intriguing connection between outflows today and in the 
past. Hence our momentum-driven wind model provides a 
concrete (albeit phenomenological) model for how outflows 
are related to the properties of star-forming galaxies across 
cosmic time. 

Our primary result in this paper is that our momentum- 
driven wind simulation is able to reproduce the observed 
kinematics of PLAs. Without winds, our simulated PLA 
kinematics do not show nearly enough wide-separation 
systems and have an average velocity width that is too 
small, in agreement with pr evious studies. A constant wind 
model, as i mplemented by ISpringel fc Hernquist I ((2003b) 
and used by iNagamine et al. I ( 2007 ) to study PLAs, im- 
proves over the no-wind case but is still unable to statis- 
tically match PLA kinematics. The primary reason for the 
success of momentum-driven wind scalings is that it ejects 
large amounts of mass from small galaxies, at relatively low 
velocities such that the gas is not overly heated. Although 
our simulations are somewhat simplistic, using a heuristic 
wind prescription, no radiative transfer, and relatively low 
spatial resolution (~ 200 pc physical at z = 3), they are 
the first to statistically match the observed PLA kinematics 
within a cosmological context. The implication is that PLAs 
represent the extended neutral gas envelopes of early star- 
forming galaxies that are driving outflows, and furthermore 
that such outflows are effective at moving large amounts of 
neutral gas into the regions surrounding early galaxies. 

Our paper is organised as follows: In Sj2]we describe our 
simulations and methodology for computing PLA absorp- 
tion. In f[3] we present our results for PLA abundances and 
present statistical tests on the kinematics of each feedback 
model compared to the observations. We also present many 
physical properties that relate PLAs to their host galaxies 
and dark matter halos. We summarise our results in §4. 



2 SIMULATIONS 
2.1 Input physics 

We simulate 8 /i _1 Mpc periodic cubic comoving volumes 
using the cosmological tree-particle mesh-smoothed par- 
ticle hydrodynamics (SPH) code Gadget-2, with mod- 
ificati ons including radiative coo ling and star forma- 
tion ISpringel fc Hernquist I 12003 j). Our version used in 



this work (described in lOppenheimer fc Pave I l2006h als o 
includes metal- line cooling 1 Sutherland fc Popita I 19931 ). 
assumes a (spatial ly-uniform) cosmic photoionising back- 
ground taken from lHaardt fc Madaul ()200lh , and includes 
a heuristic prescription for galactic outflows driven by star 
formation that we describe in the next section. We focus on 
z = 3 outputs as that is the typical redshift for the PLAs 
with well-measured kinematics. 

We choose cosmological parameters of Q m — 0.3, J1a = 
0.7, h — 0.7, as = 0.9, and Q b = .04 and generate the ini - 
tial conditions at z = 199 using the lEisenstein fc Hu I (|l999h 
transfer function. These parameters are consistent with the 
WMAP-1 results (|Spergel et al. 1 120031 ') but are somewhat 
differ ent than the latest WMAP-7 results (|Komatsu et al. I 
2008). However, these differences are not expected to affect 
our general conclusions, as the uncertainties in modeling the 
baryonic physics probably dominate over these differences. 
We use 256 3 particles each for the gas and the dark matter, 
yielding particle masses of 4.84 x 10 M@ for the gas and 
3.15 x 1O 6 M0 for the dark matter. The equivalent Plummer 
gravitational softening length is 0.625 h~ x comoving kpc, or 
223 proper pc at z — 3. These simulations resolve star for- 
mation in halos down to virial temperature of ~ 2 x 10 4 K, 
below which ambient photoionisation is expected to provide 
significant suppression of gas accretion. We note that the 
numerical resolution in our cosmological volume is better 
than that in the individual halo simulations of HSR98. 

2.2 Outflow models 

Galactic outflows appear to be related to star formation, 
but the exact relationship is unknown. While many re- 
cent studies have focused on explicitly driving winds by 
over-pressurising th e ISM via supernova heat input (e.g. 
IStinson et al~l 120061 : ICeverino fc Klypin I |2008h , it is clear 
that such processes are sensitive to scales below any realis- 
tically achievable resolution within a cosmological volume. 
Hence SH03 took the approach of explicitly incorporating 
outflows with tunable parameters, to avoid a strong depen- 
dence on physics below the resolution scale. 

In the SH03 implementation, an "outflow model" is de- 
scribed by choosing two parameters: the mass loading factor 
rj, and outflow velocity u w i n d. The mass loading factor is de- 
fined as rj = M w /Mi,, where M w is the mass loss rate by 
winds and is the star formation rate. 

Each gas particle that is sufficiently dense to allow star 
formation has some probability to either spawn a star par- 
ticle, or to be kicked into an outflow. The ratio of those 
probabilities is given by the mass loading factor 77. If a par- 
ticle is selected to be in an outflow, its velocity is augmented 
by Dwind, in a direction given by ±vxa (resulting in a quasi- 
bipolar outflow). 

When a gas particle becomes a wind particle, its hy- 
drodynamic forces are turned off until the particle escapes 
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the star-forming region and reaches a density of one-tenth 
the critical density for star formation. The maximum time 
for hydrodynamic decoupling is 20 kpc divided by the wind 
speed, ensuring that the hydrodynamic forces will almost 
always be turned off for less than 20 kpc. At the simulation 
resolution we employ, our results are not very sensitive to 
whet her or not we decouple (but see lDalla Vecchia fe Schavel 
2008, who show that it makes a large difference at very high 
resolution), but it does allow for better resolution conver- 
gence. However, this decoupling does imply that the detailed 
distribution of metals around galaxies may be inexactly 
modeled owing to the lack of accounting for hydrodynamical 
effects. For this reason, we will prefer to use the H I distribu- 
tion directly when computing DLA kinematic s, and simply 
tie the metals directly to the H I. We note that lTescari et all 
(2009), using simulations similar in physics and resolution 
to ours, explored the difference in kinematics between us- 
ing the metals directly and using the Hi distribution with 
a fixed assumed metallicity, and found that observed kine- 
matics were better reproduced in the latter case, suggesting 
that some metal diffusion or smoothing may be required. 
We follow their more optimistic case here, with the hope 
that future simulations will allow us to more robustly and 
self-consistently model the metal distribution near galaxies. 

SH03 chose M w i n d and r\ t o be constants, based on obser- 
vations available at the time ((Martin Ill99sl :i Heckman et al. I 
2000). The wind velocity v w was derived from the supernova 
feedback energy, 



1-2 

-M w v w = x^snM*, 



(1) 



where esN is the energy deposition by Type II supernovae 
per unit of star forming mass, and \ ls the fraction of super- 
novae energy that drives the wind. SH03 set x = 1> an d used 
esiv ~ 4x 10 48 erg/M0 for a Salpeter initial mass function, 
yielding a wind velocity of i> W md = 484 km/s. They set r\ = 2 
to roughly reproduce the observed present-day mass density 
in stars. SH03 showed that this wind model produces a cos- 
mic star formation history that is in broadly agreement with 
observations, and is resolution converged. We will refer to 
this as the constant wind ("cw") model, since both r\ and 
Wwlnd are independent of galaxy size. 

Improved recent observations of local starburst out- 
flows suggest that the wind speed is not c onstant, but is 
proportional to a g alaxy's circular velo city (|Martin 1 120051 ; 
iRupke et ai~ll2005l ; IWeiner et al. 1120091 ). This is suggestive 
of momentum-driven winds, as worked out analytically by 
iMurrav. Quataert. fe Thompson] (120051 ). OD06 implement 
a momentum-driven wind model in a cosmological simula- 
tion. In this scenario, the wind speed scales as the galaxy's 
velocity dispersion <r, and since the momentum deposition 
per unit star forming mass is assumed to be constant, the 
mass loading factor is inversely proportional to a. 

For our momentum-driven wind model, we employ 



&ry/JL 

CO 



(2) 
(3) 



where fi is the galaxy luminosity divided by its Eddington 
luminosity and <to is a normalisation factor. Because our 
resolution is still too low to compute the galaxy's velocity 
dispersion, we take a = yj — ^<I> from the virial theorem, 



where $ is the local potential depth where a wind particle 
is created. T he normalisation factor op is set to 300 km /s as 
suggested bv lMurrav. Quataert. fe Thompson] (12005). Val- 
ues for Jl ar e taken from observat ions: iMartin l l|2005h found 
f L ~ 2, and IRupke et al. I |2005l ) measured f L w 1.05 - 2, 
for local starbursts. In the momentum-driven wind scenario, 
lower metallicity star formation should produce stronger 
outflows owing to greater ultraviolet (UV ) flux output; w e 
account for this using the stellar models of lSchaerer I (2003). 
Specifically, we set 



h = fm x io 



-0.0029X (log Z+9) 2 ' 5 +0. 417694 



(4) 



where /l,q is randomly chosen between 1.05 — 2.0, and Z 
is the gas particle's metallicity. We term this momentum- 
driven wind model ("vzw"). In contrast with the cw model, 
Wwind and rj depend on galaxy size. 

We note that this momentum-driven wind imple- 
mentation differs fro m our more recent works (e.g. 
lOppenheimer fe Dave I 120091) that use the galaxy 
mas s to calculate a. Howe ver, as we showed 
in lOppenheimer fc Dave I l|2008l ). this choice makes 
little difference at z = 3, and only becomes critical at 
2 < 2 when large potential wells with hot gas develop. 
Furthermore, our im plementation also differs from that in 
iTescari et al.l (2009), who use a calibrated relationship to 
the halo mass to determine a. 

Finally, for comparison we also evolve a model with no 
galactic outflows (no winds, "nw"). Note that this model still 
includes thermal supernova feedback within the context of 
the sub grid multiphase ISM model of ISpringel fc Hernquist I 
(2003a). It is a poor match to many observations, but pro- 
vides a baseline to assess the impact of outflows. The suite of 
three models employed here are a subset of those in OD06. 



2.3 Distribution of neutral gas 

Here we describe how we calculate the distribution of neu- 
tral gas in each simulation. There are two significant com- 
plications for this: firs t , in th e multi-phase ISM model of 
ISpringel fc Hernquist I (|2003al) that we employ, only some 
fraction of each particle is actually neutral. Second, since 
we do not perform detailed radiative transfer, we must ac- 
count for the effects of self-shielding that can mitigate the 
strength of the photoionising flux. 

The first issue is straightforward to handle, at least to 
a sufficient approximation. In our multi-phase model, each 
particle above a critical density for fragmentation is assumed 
to consist of a cold phase at 1000 K plus a hot phase at 
10 6 K. The cold phase dominates the mass fraction, ranging 
between 84 — 100% (for primordial gas). Since we are only 
interested in the cold phase, we assume all particles above 
the fragmentation density to have 90% of their mass in the 
cold phase at 1000 K, and ignore the hot phase. Note that 
with metal cooling included, the fragmentation scale can 
depend on particle metallicity (OD06) and the cold mass 
fraction can range slightly lower, but to the level of approx- 
imation required here this is not important. In detail, the 
cold phase is probably confined to even smaller sub-particle 
scale clumps, but as long as the overall cross-section of these 
clumps is close to unity on galactic scales (as is found obser- 
vationally when examining the Hi content of galaxy disks, 
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e.g. IZwaan et all 120051 ). then the assumption of smoothly- 
distribution gas within the particle seems reasonably valid. 

Accounting for self-shielding requires more signif- 
icant approximations. The observed HI column den- 
sity distribution shows an excess of absorbing systems 
with N(HI) ^ 2 x 10 20 when the distribution is ex- 
trapolated from the low column density region pop- 
ulated by the Lyman a fore st (|Lanzetta et al. I Il99ll ; 
St orrie- Lombardi fe Wolfe I f2000). This excess has been ex- 
plained as a self-shielding effect typically causing a sharp 
transition lay er from mostly ionised regions to mostly neu- 
tral regions dMurakami fc Ikeuchi I Il990l ; iPetitieanl 1 19921 ; 
ICorbelli. Salpeter. fc Bandeira Il200ll ). Hence self-shielding 
is a critical aspect for setting the ionization level of gas in 
DLAs. 

Owing to our lack of a full radiative line transfer treat- 
ment of Lya photons, we resort to a simpler estimate for 
self-shielding. HSR98 use a pure density criterion which as- 
sumed that all the gas with a number density above 0.01 cm 2 
is fully self-shielded, while the rest is subject to the full 
metagalactic UV flux. In reality, such a threshold depends 
on the size of the collapsed object and the local strength of 
UV background, and hence should be different at each local 
position. 

Our approach is to have two criteria to identify a gas 
particle as "self-shielded": a maximum temperature and a 
minimum density. We choose a single value for each crite- 
rion within each wind model, and use the observed abun- 
dance of DLAs to constrain this value. We will demonstrate 
that the DLA abundances are not sensitive the choice of 
temperature threshold, and that the DLA kinematics ob- 
tained using different (reasonable) thresholds do not differ 
significantly. For densities, we will employ comoving densi- 
ties of pe = (40, 80)p C rit, where p CI it is the critical density 
at z = 0, as two thresholds that span a reasonable range. 
These correspond to (1000, 2000)p(,, or a baryon density of 
approximately (0.014,0.028) cm -3 (physical). 

In Table[T]we summarise our models. The names signify 
the wind model and pg. Also shown are the DLA number 
densities per unit redshift, and the number of randomly cho- 
sen lines of sight (LOS) in our DLA sample. These LOS will 
be used to investigate the kinematics of each model. Mod- 
els in boldface, cw40 and vzw80, are the ones that broadly 
match the observed number density of DLAs, as we will 
discuss further in EI3.1I while also having an observationally- 
consistent cosmic star formation rate at z ~ 3. When we 
need to compare these wind models with the no-wind case, 
we choose nw80 as a representative, although this model 
produces too many DLAs and strongly overpredicts the cos- 
mic star formation rate. Note that in order to bring the 
no-wind model into agreement with the number density of 
DLAs, we need a very high choice of pg that is likely to be 
unphysical (in the sense that this density well exceeds the 
threshold density for our multi- phase ISM model). This fur- 
ther highlights the impact that outflows have in regulating 
the amount of neutral gas in DLAs and the cosmos. 

2.4 DLA absorption profiles 

We now identify DLAs and calculate the absorption profiles 
for low-ionisation species. We first perform self-shielding cor- 
rections to the neutral gas content of individual gas particles 



Table 1. 



Name 


P0 a 


dN/dz h 


ATlos c 




-^30km/s e 


nw40 


40 


0.534 


1995 


883 


404 


nw80 


80 


0.443 


1995 


716 


310 


cw40 


40 


0.244 


1995 


1199 


672 


vzw40 


40 


0.379 


1992 


1306 


973 


vzw80 


80 


0.265 


1988 


1453 


1175 



a Self-shielding comoving density in units of the z = 
critical density. 

b Number of DLAs per unit redshift. 

c Number of lines of sight. 

d # of DLAs with Av > 20 km/s and satisfying eq. [9] 
e # of DLAs with Av > 30 km/s and satisfying eq. [9] 

as described in the previous section. To identify DLAs, we 
make an 8000x8000 Hi map projected through the entire 
simulation volume. We randomly choose 2000 pixels with 
iVm > 10 20 ' 3 cm -2 , and compute Hi Lya spectra along lines 
of sight (LOS). We then identify the individual absorbing re- 
gion along the LOS that is a DLA. In less than 1% of the 
cases, the total Nm of any single absorber along the LOS 
is below the DLA threshold, even though the total column 
exceeds it; we discard those LOS. The actual numbers of 
identified DLAs for each model are listed in Table 1. 

To obtain DLA kinematics, we attempt to mimic the 
observational procedure of identifying low-ionisation, rela- 
tively unsaturated metal lines that can trace the gas mo- 
tions. Although we track metal enrichment in our simula- 
tions, we choose not to attempt to make spectra from the 
metals directly. This is because (1) we don't accurately track 
the ionisation field, and (2) as explained above, the tracking 
of metals in outflows near galaxies is hampered by limited 
numerical resolution and explicit hydrodynamic decoupling. 
Instead, we tie the metal abundance to the Hi abundance 
with values representative of fully neutral gas. 

We choose Sill (1808A) as our representative low- 
ionisation metal tracer, which was used in HSR98 and is 
generally a good tracer for DLA kinematics. We assume a 
uniform silicon abu ndance of [Si/H] = —1, a nd ta ke the oscil- 
lator strength from lTripp, Lu. fc Savage! (| 19961 ). Since the 
ionisation potentials of Si I and Sill are 8.15 eV and 16.3 eV, 
respectively, which brackets H I, we assume that Sill is domi- 
nant and exists only in self-shielded regions and in cold ISM 
gas. This gives an abundance of nsai/^Hi = 3.24 x 10~ 6 . 
We multiply the H I optical depths along each spectrum by 
this value, along with the oscillator strength difference, to 
obtain the Sill optical depths, thereby generating spectra 
from which we can examine DLA kinematics. Our results 
are not sensitive to the specific abundance and ionization 
choice, because we are primarily interested in the redshift- 
space extent of metal systems rather than the strengths of 
the systems themselves. 

Example randomly-chosen spectra are shown in Fig- 
ure [T] Each panel contains 5 LOS; the locations (in our 
8000 x 8000 grid) and Hi column densities are indicated 
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Figure 1. Example low-ionisation metal line (i.e. Sill) spectra 
drawn from our simulations. The metal optical depths are tied 
to Hi as described in the text. Each colour shows a different 
LOS; solid and dotted lines show smoothed (19 km/s tophat) 
and unsmoothed versions, respectively. 

along the right side. Such spectra are analyzed as described 
in ^3.21 in order to obtain kinematic statistics for each DLA. 

2.5 Identification of galaxies and dark matter 
halos 

Besides DLA kinematics, we are also interested in studying 
the relationship between DLAs and galaxies in our simula- 
tions. Hence, we must ide ntify galaxies in ou r simulations 
and we do so using SKIEQ (|Katz et al. Ill996al ). which iden- 
tifies density watersheds using spline kernel interpolation, 
and then finds surrounding groups of gas and star particles 
that are gravitationally self-bound within these density wa- 
tersheds. We apply SKID to all the star particles and gas 
particles that satisfy Pgas//0gas > 10 3 and T < 30, 000K, 
and all particles (regardless of temperature) above the mul- 
tiphase threshold of « 1000p S a S . 

Figure [2] shows the galaxy baryonic mass function for 
our three simulations We take our g alaxy mass resolutio n 
limit to be 32m S p H = 1.55 x 10 7 M Q (jFinlator et al. Il2006t) , 
shown as the vertical line. This limit appears to be conser- 
vative, as the mass function doesn't turn over severely un- 
til much lower masses, but detailed resolution testing shows 
that the individual galaxy star formation histories are poorly 
converged below this mass owing to stochastic effects and 
the fact that the densities cannot be adequately resolved. 
The inflection point in the no-wi nd case we bel ieve corre- 
sponds to the filtering mass (e.g. iGnedinl feoOO). which is 
the mass around which galaxies are significantly suppressed 
by metagalactic photoionization. Galaxies below this mass 
formed prior to turning on our spatiall y-uniform ionizing 
back ground instantaneously at z = 9 (|Haardt fe Madaul 
l200ll ). A more realistic radiative transfer simulation would 
not yield such instantaneous reionization, and would likely 
not produce such a feature. The wind models have much 
lower stellar masses for a given halo mass, and hence the 

1 Spline Kernel Interpolative DENMAX; 

http://www-hpcc.astro.washington.edu/ tools/skid. html 



100.00 




0.Q1 1 i \ i i . ... . \ 

10 7 10 s 10 9 10 

Log galaxy mass (M@) 



Figure 2. The mass function of gravitationally bound groups 
found by SKID for the three models: no winds (nw, red solid), 
constant winds (cw, green dashed) and momentum-driven winds 
(vzw, blue dotted). We resolve galaxies down to a baryonic mass 
of 1.55 X 10 7 Mq, indicated by the vertical dotted line. 

effects of filtering are relegated to much smaller masses, at 
or near our resolution limit. While these effects on the mass 
function are interesting, they are not immediately relevant 
to DLAs, so we do not discuss them further. 

To identify dark m atter halos we use a friends-of- 
friends(FOF) algorithm (|Davis et al. I Il985l ). which finds 
groups of dark matter particles by linking neighbouring par- 
ticles within a given linking length. We choose the linking 
length to be the interparticle distance at one-third of the 
virial overdensity p V i r /p, which is the local density at the 
virial radius for an i sothe rmal sphere, and we take p V i r from 
iKitavama fc Sutol l|l996l ). To define the virial mass, virial 
radius, and the circular velocity of the dark matter halo we 
refine the group using a spherical overdensity (SO) criterion. 
In SO, the halo centre is set to be the location of the most 
bound FOF particle. Then we expand the radius around this 
centre until the mean overdensity inside the radius equals 
Pvir/p. We define this radius to be R v i T , the mass within R v i T 
to be M v ir, and the circular velocity v c — \J GM v - a / R v - a . 

In the end, we obtain a sample of about 2000 DLAs and 
~ 3000 galaxies within each of our simulations. We will now 
analyze these systems to understand the impact of outflows 
on DLA kinematics and the relationship between DLAs and 
galaxies. 



3 PHYSICAL PROPERTIES OF DLAS 
3.1 DLA number densities 

The observed redshift-sp ace abundance of DLAs, dN/dz = 
0.26 ± 0.05 (z = 3 ; IStorrie-Lombardi fc Wolfe! l2000l ; 
IProchaska et al. II2005I) . provides a key constraint for DLA 
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models. In our case, we use it to constrain our self-shielding 
criteria, pg and Tg. Figure[3]shows our derived abundance for 
each wind model assuming different pg and Tg. Remember 
that we assume that all gas with a density greater than pg 
and a temperature less than Tg is fully self-shielded. Overall, 
changing Tg has little effect on the abundance, while pg has 
a large effect. 

The no-wind (nw) model has too much neutral gas to 
match the observations within a large range of choices for 
pg. For the two wind models, we choose two self-shielding 
criteria, pg — 40 and pg = 80 for the cw model and the 
vzw model, respectively, to approximately match the ob- 
served dN/dz. A density of pg = 40 corresponds to lOOOpf, 
and riH = 0.014 cm -3 in primordial composition gas, which 
was adopted in HSR98 as their shielding criteria. So the cw 
model has the same density threshold value as HSR98 while 
the vzw model has one twice as large. Our resolution and 
inclusion of metal line cooling likely affects the choice of 
threshold values. 

Figure [4] shows HI column density maps projected along 
the x-axis for the nw80, cw40, and vzw80 models. We see 
substantial differences in the neutral gas distributions pro- 
duced by these three different feedback implementations. 
Without winds, the DLA-absorbing gas is highly concen- 
trated within star-forming regions. In the constant wind 
(cw) case, the absorption is somewhat more extended, but 
the high outflow velocities generally cause the outflowing gas 
to be heated, thereby lowering their neutral fractions and 
hence their densities. The momentum-driven wind model 
(vzw) produces the greatest extent of high-column density 
gas. This foreshadows our result that the vzw model will 
yield the largest DLA kinematic widths. 

The more extended gas distribution causes the vzw 
model to be most sensitive to pg. This subtle dependence 
arises because the abundance is sensitive to cross-section, 
since an extended distribution results in more gas lying near 
the threshold density. So the curves in Figure [3] partly re- 
flect the different gas topologies that result from the different 
outflow models. In contrast, the sensitivity to temperature 
threshold is negligible because most of the gas for any rea- 
sonable choice of pg lies at ~ 10 4 K, below which we truncate 
radiative cooling in our simulations. 

As an aside, we note that the filamentary gas struc- 
tures that are responsible for feeding galaxie s with fresh 
fuel via the cold mode accretion scenario (e.g. iKeres et al.l 
120051 ; iDekel et ai1l2009f ) are generally not sufficiently dense 
to produce DLA absorption, and instead typically have 



Only in the vicinity of galaxies does the 



gas starts to self-shield, and the DLAs tend to be confined 
to even denser gas within and around galaxies. 



3.2 DLA kinematics 

We now examine the primary target of our investigation: 
The kinematics of simulated DLAs and how they compare to 
observed kinematics. PW97 developed a suite of statistics to 
descibe the one-dimensional distribution of metal absorption 
lines within a single DLA. To generate these statistics from 
our simulated DLAs, we adopt the following procedure: 

(i) We select DLAs and then generate Sin absorption 
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Figure 3. DLA abundances for each model assuming various self- 
shielding criteria. The temperature dependence is weak, while the 
density dependence is strong. The red arrows indicate the value 
we choose for pg. Note that the values in this figure are slightly 
larger than the values in Table [T] because we used 2000x2000 
maps for this figure and 8000x8000 maps for our main results as 
presented in Table [T] 



lines by tying Si 11 to H I as described in £12.41 

(ii) We smooth the absorption lines with a 19 km/s 
tophat window, which is equivalent to the 9-pixel tophat 
smoothing done in PW97. 

(iii) From this profile, we calculate four key quantities, 
as defined in PW97: The velocities of the largest (i> p k) and 
second largest peak absorption (i>2 P k), the median absorber 
velocity (iw,d)> and the system velocity width (Aw) defined 
by 



«9S - V5\ 



(•>) 



where 1)95, vs, and u mc d represent the positions at which the 
optical depths are 95%, 5%, and 50% of the total optical 
depth. We note that we compute these quantities directly 
from the Sin optical depths, while PW97 compute them 
from apparent optical depths obtained from the flux; we 
expect these quantities to be similar, particularly for the 
typically unsatured metal lines in DLAs. 

(iv) From these quantities, we calculate the DLA kine- 
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Figure 4. The column density of neutral hydrogen (HI maps) for 
the no wind model (nw80, top panel), the constant wind model 
(cw40, middle panel) and the momentum driven wind model 
(nw80, bottom panel). Every pixel over 2x 10 20 cm — 2 , i.e. damped 
absorption, is plotted in red. All lengths are comoving. 



matic statistics from PW97: 

^mcd fmean | 



fm 

/cdg — 
/2pk = ±' 



(Av/2) 

|^pk — V mean | 
|f2pk — Vmean| 

(Aw/2) 



(6) 
(7) 
(8) 



where the plus sign for the two peak fraction, /2 P k, holds if 
the second peak is between the mean velocity and the first 
peak velocity; otherwise the minus sign holds. If there is no 
second peak, we take the edge- leading fraction, / e d g , for the 
second peak fraction / 2pk . Here, « mC an = |(^9S + vs). 

To avoid saturation effects that blur the kinematic infor- 
mation, and to exclude noise contamination in the obser- 
vations, we again follow PW97 by only using profiles with 
peak intensities / p k in the range 



0.1 < < 0.6 
la 



(9) 



where Jo is the continuum around the absorption line, or 
peak optical depths between 0.5 and 2.3. This removes a 
significant portion of the total DLA sample, roughly 60% 
for nw, 45% for cw, and 30% for vzw. The variance amongst 
models in the acceptance fraction based on this criterion 
suggests that this may be another way to constrain outflow 
models. 

The most crucial statistic is Aw, the system velocity 
width. It represents the velocity-space extent of the dense 
neutral absorbing gas, and hence encodes information about 
internal motions within the ISM as well as any inflow or 
outflow- induced motions. The other statistics, / mm , /cdg, 
and /2 P k, turn out to be less discriminatory, but we include 
them for completeness. / mm measures the skewness of the 
overall absorption in the DLA; a symmetric distribution of 
optical depths would yield / mm = 0. The edge-leading test 
/ c dg would be if the strongest absorption is at the kine- 
matic center, but is large if the kinematics are dominated by 
rotation or infall where the strongest absorption occurs at 
large velocities from the center. The 2-peak test /2 P k is de- 
signed to distinguish between rotation and infall: in the case 
of rotation, the second peak is expected to be on the same 
side as the first (and hence yield a positive value), whereas 
spherical and symmetric accretion would produce the peaks 
on opposite sides, yielding a negative value (PW97). While 
these statistics were devised to distinguish between simple 
scenarios, the complex interplay between infall, outflow, and 
rotation within a fully hierarchical context precludes such 
straightforward interpretations. Hence we focus on the dis- 
tributions of these statistics among DLA samples (both ob- 
served and simulated), and use Kolmogorov-Smirnoff (K-S) 
tests to characterize their (dis)agreement. 

In Table [1] we summarise the number of simulated ab- 
sorption lines, A^ampic, taken from each model. For the ob- 
serv ations, we will compare to 46 observed lines of sight 
from lProchaska fc Wolfel |200ll) . Note that their paper only 
presented the distribution of Aw; the other quantities were 
kindly provided by X. Prochaska (private communication). 

Before we compare the simulations with the observa- 
tions, we must consider the issue of velocity resolution ef- 
fects in the observed DLA sample. The 9 pixel smoothing of 
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PW97 effectively provides a minimum of Av m in ~ 20 km/s, 
and indeed all of the observed DLAs have Av > 20 km/s. 
But HSR98 adopted a minimum threshold of Av m in = 
30 km/s to avoid any possible incompleteness of the DLA 
sample in the range of Av = 20 — 30 km/s. Unfortunately, 
this threshold can bias the statistics, as it turns out that 
there are a significant number of smaller halos that host 
DLAs having Av within this range. To test this incomplete- 
ness issue, we present results for both velocity thresholds of 
20 km/s and 30 km/s. However, we note that the observed 
sample may be significantly incomplete for the 20 km/s 
threshold, and may not be a fair comparison to the mod- 
els, hence we prefer to compare to the 30 km/s case. As it 
happens, our main conclusions are unaffected by this choice. 

Figure[5]plots histograms of the four PW97 statistics for 
our three outflow models, along with their K-S test values as 
compared to the observed statistics, with Af m m = 30 km/s. 
The upper left panel demonstrates the central result of our 
paper: Models with outflows provide a much better match 
to PW97's velocity width statistic than our no- wind case. In 
particular, the momentum-driven wind simulation provides 
an excellent match to the data with a K-S acceptance level of 
0.39, although the constant wind case cannot be definitively 
ruled out. The rest of the paper will mostly be devoted to 
understanding the physical origin of the differences among 
model for the velocity width statistic Av. We will see that 
the variations arise because the differences between wind 
models are most significant for small halos. 

The other PW97 tests generally show acceptable K-S 
values for all models when compared to the data, at least 
for Av m i n — 30 km/s; at An m i n = 20 km/s again the vzw 
case is clearly favored as seen in Table [2] The median-mean 
statistic / mm (Figure O top right) shows that the cw model 
produces slightly more skewed DLAs than the vzw and nw 
cases, resulting in marginally poorer agreement in the K- 
S test, but not so significant as to reject the model. The 
edge-leading statistic /edge (bottom left) tells a similar story, 
that nw produces the fewest edge-leading profiles and cw the 
most, but none of the models can be excluded by the data. 
The two-peak test (bottom right) shows that most galax- 
ies show rotation signature in their center and interception 
with infalling or outflowing gas is rare, but it is more fre- 
quent for the wind models. This last statistics shows the 
most disagreement between the simulations and the data, 
and none of the models are obviously favored. It is unclear 
why this discrepancy is strongest, since it is a difficult statis- 
tic to interpret. However, we note that the thick disk model 
favoured by PW97 shows a similar distribution, and hence 
is expected to have roughly the same K-S value. 

Our no-wind results are significantly different than 
HSR98, who found broad agreement with their simulations 
(without strong outflows) and the data. The discrepancy 
arises because they extrapolate the relationship between 
DLA cross-section and halo mass to low halo masses in or- 
der to obtain cosmologically-based statistics for their DLA 
velocity widths. In ^3.4.31 we calculate this relationship di- 
rectly from our simulations, and show that their extrapola- 
tion is not valid. Specifically, our no-wind simulation shows 
considerably increased cross-section over such an extrapola- 
tion. This is the main effect that causes our no-wind case 
to be in poor agreement with data, as compared to HSR98. 
We note that our numerical resolution is similar to, in fact 
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slightly better than, HSR98. Since HSR98 simulated indi- 
vidual galaxies while we employ a full cosmological simula- 
tion, we are able select DLAs based upon their absorption 
cross-sections to match the observed abundance of DLAs 
(which HSR98 could not do). This enables us to more accu- 
rately characterize the simulated DLA population and more 
robustly compare it to observations. 

HSR98 also noted that their results depended on the 
assumed self-shielding criterion. In Table [2] we present the 
results for different self-shielding criteria as well as for dif- 
ferent minimum velocity thresholds. The K-S tests indicate 
that the vzw model is consistent with the observations for 
either choice of the threshold density, while the no-wind 
case is strongly discrepant in Av regardless of threshold 
choice. Hence reasonable variations in the self-shielding cri- 
teria do not significantly impact our overall conclusions. 
Hence while it certainly possible to perform more sophis- 
ticated s elf-shielding models e ven without radiative trans- 
fer (e.g. iPopping et al. 1 |2009| L_ or even do the fu ll post- 
processing radiative transfer i Pontzen et al. ||2008| ). we do 
not expect that this will strongly impact our results. 

In summary, galactic outflows appear to be capable of 
reconciling CDM-based models of galaxy formation with ob- 
served DLA kinematics. This is the most important conclu- 
sion from our work. Such statistics also provide constraints 
on wind models, and in the limited tests conducted here we 
favour momentum-driven wind scalings over constant wind 
scalings. We now investigate these trends more deeply by 
studying the properties of galaxies and halos that give rise 
to DLAs, in order to understand the physical reasons behind 
the success of outflows in general, and the momentum-driven 
wind scalings in particular. 

3.3 Distance to DLA host galaxies 

In our simulations, DLAs are generally associated with 
galaxies. This is believed to be true in the real Universe, 
although it is usually only possible to image the host galax- 
ies at low redshifts, where galaxies can be separated from the 
bright background quasar. On the other hand, DLAs seem 
to deviate from well-established pro perties of galaxies suc h 
as the Kenni cutt-Schmidt relati on (|Wolfe fc Chen I l2006h . 
Furthermore, IWolfe et al. I l|2008l ) identified a bimodality in 
DLAs, in which they divide into systems with high and low 
[CII] 158 fim cooling rates, with the former having proper- 
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Figure 5. Distributions of the four statistics from PW97 (shaded histogram) compared with the nw80 (red), cw40 (green), and vzw80 
(blue) models. A minimum velocity threshold of 30 km/s is assumed. Numbers in parenthesis are K-S test values comparing each model 
with observations, summarised in Table[2] The no- wind simulations is strongly disfavored, and statistically the vzw model provides better 
agreement than cw. 



ties similar to Lyman Break Galaxies and the latter perhaps 
arising in a different sort of population. Hence the relation 
between DLAs and galaxies remains uncertain. 

Here we examine some basic properties of the relation- 
ship between DLAs and galaxies in our simulations. To con- 
nect DLAs to their host galaxies and their host dark mat- 
ter halos, we need to define the line-of-sight positions of 
DLAs. Since it is somewhat arbitrary to define the position 
of diffuse gas, we consider two definitions: x$o = the line- 
of-sight position where 50% of the optical depth is reached; 
and a; a ve = 0.5(2:30 + X70), where £30 and £70 are defined 
analogously to x$o. We also define an uncertainty on the 
line-of-sight position as Aa; = \x 5 o — x mc \. 

In the vast majority of cases, Aa; is quite small, so the 
distance is fairly unambiguous. The few percent of cases 
where it is not are typically caused by two or more separate 
gas clumps along the line of sight contributing to a single 
system. For the remainder of this paper, we will exclude 
cases where Aa; > 0.01 in units of the simulation box length 
(i.e. 28.7 proper kpc), since these systems are difficult to 
relate unambiguously to galaxies. Specifically, this excludes 



1.5%, 6.9%, and 3.3% of DLAs in the nw, cw, and vzw cases, 
respectively. We also examined a tolerance criterion of Aa; = 
0.001, with negligible difference in the final results. 

The top panel of Figure [6] shows the projected dis- 
tances, -Dproj, of DLAs versus their column densities for 
all three models. In general, the median distance is around 
~ 2 kpc (distance units are physical unless noted) in all 
models. There is a weak trend for higher-column systems to 
arise closer to galaxies. Almost all systems in all cases arise 
within 10 kpc of a galaxy. Hence in our simulations, DLAs 
generally arise in the extended neutral gas that are present 
in and around galaxies at these epochs. 

The middle panel shows the cumulative distribution of 
projected distances, P(> D plo j), for each wind model. The 
no-wind case shows a qualitatively different behavior than 
the wind runs: It shows somewhat larger typical distances 
at small separations, but beyond > 3 kpc the nw case shows 
very few DLAs, while the wind models possess a significant 
tail to high projected separations. Hence about 95% of the 
DLAs in the no wind case have distances smaller than 5 kpc, 
but outflows substantially puff out the neutral gas as shown 
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Figure 6. The projected distance, -D pro j, from the host galaxy 
(top panel), its cumulative distribution (middle panel), and 
the cumulative distribution of D pro] /D, which characterises 
the angular distribution (bottom panel) for the all three mod- 
els as labelled. The purple dotted curve has the form 1 — 
y/l — (Z) pro j/£)) 2 , which is what one would expected for an 
isotropic distribution of small blobs. The no- wind case is clos- 
est to having a central blob topology, while winds (particularly 
vzw) moves the gas distribution closer to isotropic. 



in Figure 13] so that in the vzw case 95% of DLAs are within 
10 kpc. Since the mass loading factor is larger for small halos 
in the vzw model than in the cw model, the momentum- 
driven wind model is most efficient at puffing out the neutral 
gas from the more numerous and gas-rich small systems. 

One can quantify the topology of DLA absorption by 
examining the ratio of projection distance D pro j to real dis- 
tance D, shown in the bottom panel. First, let us consider 
some simple illustrative cases: If DLAs come from small 
blobs that are distributed isotropically from host galaxies, 
the cumulative distribution should follow the purple dotted 
line, P(< D ploj /D) = 1 - y/l - (D ploj /D) 2 . This is likely 
similar to a scenario in which DLAs arise from randomly- 
dis tributed smal l disks within a galactic halo, as forwarded 
bv lMaller et al. I l|200ll ). Conversely, if all DLAs are located 
in the cross-section plane, then the cumulative distribution 



should be a step function at a value set by the angle of the 
plane relative to the line of sight. A special case of this is if 
the DLA owes to a single, spherical gas blob centered on a 
galaxy ("central blob"), in which case the distribution will 



be a step function around unity since D p 



D. 



Now let us examine the actual distributions. The three 
models generally lie between isotropically-distributed blob 
and central blob geometries, not surprisingly since these 
are the extreme cases. In more detail, the no-wind model 
lies distinctly more towards central blob, because most of 
the DLA cross-section is concentrated in the central galaxy. 
Conversely, the wind models follow more closely to the 
isotropically-distributed blob case, with the vzw case more 
so than the cw. This means that outflows change the sim- 
ple centrally-dominated neutral gas topology of the no- wind 
model to a more complex and extended gas topology. This 
is consistent with the general impression one gets from ex- 
amining the H I images in Figure [4] 

These topological trends directly correlate with the 
DLA kinematics examined in i]3.2l The more the winds push 
out gas, the larger the velocity widths. The extra velocity 
can arise owing to the kinematics of the outflowing gas it- 
self, or else the fact that the gas velocity field typically differs 
more from the galaxy systemic velocity as one goes farther 
out into the halo; we dissect these scenarios further in the 
next section. In either case, the effectiveness by which the 
momentum-driven winds push out gas into the halo, owing 
largely to the high mass loading factors in small galaxies, is 
critical for yielding its agreement with DLA kinematics. 



3.4 Relations between DLAs, host galaxies, and 
host dark matter halos 

3.4-1 The environment producing high-Av DLAs 

As pushing out the gas distribution into the surrounding 
halo seems to be an important aspect of reproducing DLA 
kinematics, we now investigate how halo environments im- 
pact the kinematics of DLAs. Figure [7] shows the relations 
between the velocity widths of DLAs, Av and the mass of 
nearest galaxy M ga i (top panel), the distances from the near- 
est galaxy D (second panel), the circular velocity of the 
host dark matter halo v c (third panel), and the mass of 
the dark matter halo Mhaio (bottom panels) for the no-wind 
(red points) and the momentum-driven wind (blue points) 
models. The constant wind case is not shown for clarity, as 
its trends are intermediate between the two displayed cases. 

The Mg a i — Av plot shows that Av does not have a sim- 
ple relation to the mass of the host galaxy. The nw model 
shows no significant trend, while the vzw model shows an in- 
teresting bimodal distribution for high velocity width DLAs, 
in which the nearest galaxies tend to be the most massive 
or the least massive. We will see below that these large-Au 
systems tend to occur in massive halos, and hence we can 
deduce that the high-velocity tail in the vzw model occurs 
primarily around either massive (likely central) galaxies or 
low-mass galaxies within a massive halo. Overall, more of 
the absorption seems to come from low-mass galaxies in the 
vzw case, showing that the high mass loading factor in small 
systems produces significantly more DLA-absorbing cross 
section from these galaxies. 

The D ~ Av plot shows that the highest velocity widths 



© 0000 RAS, MNRAS 000, 000-000 



12 Hong et al. 



"°2 
a 



S03 



M 1 1 1 1 l 1 1 1 1 I 
No wjrjds 




2 2 

2. « 



^- o 

t» o 

E 

w o 

> — 



-| 1 1 1 1 1 1 I- 



H — 1 — 1 — 1 — I — 1 — t 1 — i- 



~~ i — 1 — 1 — 1 — r 
M-D winds 



-| — i — i — i — | — i — i — t— 




H — 1 — 1 — 1 — I — 1 — t 1 — H 




200 
Av (km/s) 



400 



200 
Av (km/s) 



400 



Figure 7. The velocity width of a DLA, Av, versus the masses 
of the nearest galaxy (top panel), the distances from the nearest 
galaxy (second), and the mass (third) and circular velocity (bot- 
tom) of the host dark matter halo. Red points show the nw80 
model, blue points the vzw80 model. Histograms along the j/-axis 
show the distribution of DLAs along that axis, with a tick mark 
indicating the median value. In the bottom panel, the dotted line 
shows Ati = v c ; DLAs to the right of this have "supergravitational 
motion" , which is more common in vzw compared to nw. 



occur relatively far out from the host galaxy, typically 
2 — 10 kpc away. The overall median distance is approxi- 
mately 2 kpc for the vzw case and 3 kpc for nw (shown 
as the tick marks above the vertical histograms), but most 
of the DLAs with Aw > 150 km/s lie above these values. 
This means that if one wants to obtain large Av values, one 
has to arrange a significant amount of DLA cross-section at 
relatively large distances from the host galaxy. At smaller 
distances, the internal kinematics of galaxies is insufficient 
to produce large velocity widths. The vzw model, relative 
to the nw case, clearly produces more DLA-absorbing gas 
distributed farther away from galaxies. 

The Mhaio — Aii plot (third panel) shows an envelope 
of increasing Av in larger halos. In the no-wind case, the 
envelope occurs at roughly an equality between v c and Av, 
showing that with only gravitationally-induced motion, it is 
rare to produce kinematics larger than the halo circular ve- 
locity. Conversely, the vzw case shows significant numbers 
of DLAs with widths greater than v c . This is of course key 
to producing large velocity widths, and we will explore this 



further in the next section. Note that a massive halo is a nec- 
essary but not sufficient condition for large velocity widths, 
as massive halos can host DLAs with small Av. These trends 
are mirrored in the circular velocity panel, which encodes in- 
formation about the halo radius but otherwise shows very 
similar trends. 

From the above analysis, along with visual impressions 
from Figure [4j we can assemble a scenario for DLA ab- 
sorption: Small velocity width DLAs can be produced in 
all halos, but large velocity widths require a massive halo 
environment. In the momentum-driven wind scenario, this 
environment is one where a massive central galaxy is sur- 
rounded by numerous small galaxies in a clustered envi- 
ronment. The DLA then occurs at intermediate distances 
away from the nearest galaxy, which can often be the largest 
galaxy or one of the smaller systems. The wide DLA veloc- 
ity widths are then produced either by sightlines intersect- 
ing multiple gal axies within a dense region of protogalactic 
clumps (HSR98; lMaller et al. Il200ll ). or by dense gas puffed 
up by outflows. The former can occur without winds, but 
the latter requires momentum-driven outflows that produce 
a puffier neutral gas distribution and raise the Hi cross- 
section of small galaxies within larger halos, thereby yielding 
the correct fraction of wide systems. 



3.4.2 Supergravitational motion 

In the no-wind model there is no kinetic energy injection 
into the neutral gas beyond that provided by gravity. High 
velocity width DLAs with Av > v c can still arise owing to 
merging or infalling motion of neutral gas, as in HSR98, but 
this is rare. Conversely, outflows provide substantial non- 
gravitational kinetic energy, which is reflected in the signifi- 
cantly larger fraction of DLAs that have Av > v c as seen in 
Figure [7] (fourth panel). We call the kinematics induced by 
non-gravitational energy "supergravitational motion". The 
amount of supergravitational motion (Av > v c ) therefore 
quantifies the contribution of wind feedback to the kine- 
matics of the neutral gas. 

Figure [8] shows the probability distribution p(Av/v c ) 
for each model. The probability distr ibution for the no- win d 
case is consistent with previous work (|Haehnelt et al. 2000). 
The peak of Av/v c is around 0.5 and the fraction of systems 
with supergravitational motion is only 5%. The no- wind case 
shows a strong maximum at v/v c sa 0.5, and rapidly drops 
off to higher v/v c . The wind models show a qualitatively dif- 
ferent behaviour, with a less obvious maximum and a longer 
tail to high v/v c . The fraction with supergravitational mo- 
tion for the cw model and the vzw model are 13% and 22%, 
respectively. Without such high fractions of supergravita- 
tional motion, the models cannot reproduce the incidence 
of wide DLA velocity widths. Hence the distribution of gas 
by outflows, together with small-scale galaxy clustering, are 
critical for explaining DLA kinematics. 



3.4-3 DLA cross-section of dark matter halos 

To quantify the increase in the puffiness of the gas distri- 
bution in the wind models that is responsible for the wider 
velocity widths, we calculate the DLA cross-sections for each 
dark matter halo. 
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Figure 9. The DLA cross-section a (top panels) and the ratio fj, = a/TrR^ ir (bottom) as a function of circular velocity v c for the three 
wind models. Top axis is labeled by th e approximate halo ma ss. Dashed lines in top panels indicate power-law fits to v — c > 80 km/s 
cross sections. nw80 is consistent with iGardner et al. I dl997a ") for large halos, and cw40 is consistent with Nagamine et al. (2004); in 
both cases, there is extra absorption at low v c relative to an extrapolation from high-D c . vzw80, in contrast, shows a pure power law 
with index = 3.47 down to the smallest halos. The other two models have a larger cross-section for small halos, which is one reason 
why they produce too many 1ow-Ad DLAs. 
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Figure 8. The probability distribution p(Av/v c ) for each 
wind model. The fraction showing supergravitational motion, 
i.e.Av/v c > 1, are 5%, f3%, and 22% for the nw model, the 
cw model, and the vzw model, respectively. 

The upper panels of Figure [9] show DLA cross-sections 
a in comoving kpc 2 for the nw80, cw40, and vzw80 models. 
We calculate the cross-section from the projected neutral 
hydrogen maps, including only those pixels that produced 



DLAs with Av > 20 km/s. Each DLA produces 3 points 
on this plot, for projections along each of the cardinal axes. 
The relation between a and v c has traditionally been char- 
acterized as a power law with a slope we call fi, though as 
is evident from the plot a pure power-law is only a good de- 
scription for the vzw80 model (upper right). We note that if 
the cross-section w as proportional to halo ma ss, then /3 = 3, 
since M h oc v\ (e.g. lMo. Mao, fc Whltelll998h : this turns out 
to be very roughly the case, though the departures from this 
are key for understanding DLA kinematics. 

The no-wind case shows a relatively shallow power law 
(i.e. low /3) that results in large contributions to DLA cross- 
section from small v c systems. This causes a mean velocity 
width that is too small when the overall DLA abundance 
is mat ched. For comp a rison, w e show the fit from the mod- 
els of IGardner et al. I (|l997al ). which also did not include 
winds but had much lower spatial resolution. Gardner et 
al. matched the observed DLA abundance, but in fact the 
cross-section in small systems was underpredicted compared 
to our current nw model. Hence, using higher resolution sim- 
ulations without winds tend to exacerbate the disagreement 
with observed DLA kinematics. 

For the constant wind model, the power law index is 
/3 — 3.88 fo r v c > 80 km/s. This w ind model is identical to 
that used in lNagamine et al~l l|2004l ). but our slope is slightly 
larger than theirs because we fit the slope for v c > 80 km/s 
only. We choose this cutoff to highlight where the power 
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Figure 10. The cumulative global DLA cross-section (top) and 
the differential global DLA cross-section (bottom). 



law breaks in the cw model, especially in relation to the 
vzw case. For smaller values of v c , the behaviour of the cw 
model tends towards that of the no-wind case. Hence the 
constant-wind case is able to puff out gas to some extent, 
but still produces significant absorption in low-v c systems 
that causes it to match the observed kinematics less well. 

Unlike the other two models, the vzw model shows 
no break in v c vs. cross section, with an overall slope of 
/3 — 3.47. This low cross-section at small v c is critical for 
repr oducing the observed kinem atics. As argued in PW97 
and iProchaska fc Wolfel (|200ll ). most CDM galaxy forma- 
tion models fail to reproduce DLA kinematics because they 
predict too many small baryonic structures. The key aspect 
of the vzw model is that it has a high mass loading factor r\ in 
small galaxies, which lowers the cross-section of the smaller 
systems by ejecting more of its material. Indeed, the reduc- 
tion in cross section is roughly equal to (1 + rj) in smaller 
halos as expected from simple gas-removal arguments: cw is 
lower than nw by ~ x3, and vzw is increasingly suppressed 
to smaller systems since rj oc v~ x . 

The lower panels of Figure [9] show fi = a/(TvR^ iT ), i.e. 
the ratio of the DLA absorption cross-section to the pro- 
jected area of the dark halo. For the no wind case ji is al- 
most constant and cuts off rapidly for the smallest halos, 
and the cw case shows qualitatively similar behavior, low- 
ered by ~ x3 and with somewhat more scatter. In contrast, 
in the vzw wind model /i monotonically increases with halo 
circular velocity, and hence there is relatively more DLA 
cross-section at larger circular velocities that can give rise 
to high Av values. 

In Figure [10] we examine the global cross-sections to 
quantify which halo masses are predominantly responsible 
for DLA absorption, plotted cumulatively (top) and differen- 
tially (bottom) versus halo mass. We define the cumulative 
global cross section to be 



The cumulative global cross-section reiterates the re- 
sult that the vzw model is very efficient at suppressing DLA 
absorption in small halos relative to larger ones. The cw 
model produces comparably high cross-sections in large ha- 
los to vzw, but it does not suppress the low-mass halo cross- 
section, and instead its suppression occurs more at interme- 
diate masses. 

In summary, the combination of both efficiently sup- 
pressing absorption at low masses and having high ab- 
sorption at large masses relative to the no-wind case is 
what makes the momentum-driven wind model successful 
at matching the DLA kinematics. The constant wind case 
produces the latter requirement but less so the former, mak- 
ing it intermediate between the vzw and no-wind cases. 



4 CONCLUSIONS 

We examine the kinematic and physical properties of 
DLAs using cosmological hydrodynamic simulations includ- 
ing three different heuristic prescriptions for how galactic 
outflow properties are related to their host galaxies: A model 
without outflows (nw), a constant wind (cw) model where 
the wind speed of 484 km/s and mass loading factor of 
2 are independent of galaxy properties, and a model with 
momentum-driven wind scalings (vzw) where the wind speed 
scales as the circular velocity and the mass loading scales in- 
versely with it. Our no-wind model shows results consistent 
with previou s studies of hierarchic al models without strong 
outflows fe.g. lPontzen et al. Il2008h . and dramatically fails to 
explain the observed kinematics of DLAs, particularly their 
systemic velocity widths as traced by low-ionization metal 
lines. This contrasts with HSR98 who concluded that merg- 
ing of protogalactic clumps (without winds) can explain the 
kinematics of DLAs; in our case, the high frequency of wide 
lines cannot be reproduced when such a model is placed in a 
cosmological context. Hence without non-gravitational mo- 
tion, it appears that DLA kinematics are not reproducible 
within modern hierarchical structure formation models. 

We then explored DLA kinematics including outflows. 
The central result of this paper is that our momentum- 
driven wind model provides a very good match to all ob- 
served DLA kinematic measures, with the possible exception 
of the 2-peak test that no model matches well. Meanwhile, 
a model with constant wind spe eds and mass loading fac- 
tors (as in lNagamine et al. 1120071 ') is only marginally consis- 
tent with data, and a no- wind case is in strong disagreement 
with data, particularly in the distribution of DLA velocity 
widths (Av) that have historically been difficult to repro- 
duce in CDM-based models. This shows that galactic out- 
flows can reconcile DLA kinematics within current galaxy 
formation models with observations. Moreover, DLA kine- 
matics could potentially provide interesting constraints on 
the properties of galactic outflows at high redshifts. 

The abundances of DLAs may also provide interesting 
constraints. Our no-wind case could not reproduce the ob- 
served abundances without a very high and likely unphysi- 
cal threshold density for self-shielding. Meanwhile, both our 
wind models can match the abundance of DLAs with rea- 
sonable choices for this parameter. The previous study of 
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iGardner et"al~l (|l997ah with no winds matched the observed 
abundance by an incorrect extrapolation of the DLA cross- 
section versus circular velocity from low resolution simu- 
lations. Our higher-resolution simulations show that such 
an extrapolation is not valid in the the no-wind or (to a 
lesser extent) the constant wind case, since these models 
produce excess cross-section at low-« c over such an extrap- 
olation. Still, given the crudeness of our self-shielding crite- 
rion, it is difficult to assess the robustness of this result. We 
note that the s ophisticated radiative transfer simulations of 
Pon tzen et al. (2008) which include self-consistently driven 
outflows are able to match DLA number densities (as well 
as metallicities), but they did not match the kinematics. 

To understand why outflows help with DLA kinemat- 
ics, we investigated the properties of DLA absorption as a 
function of galaxy and halo properties. A consistent story 
emerges that momentum-driven winds most efficiently ex- 
pel gas into a more extended, clumpy distribution around 
galaxies, particularly for small galaxies living in dense en- 
vironments. As a result, the DLA cross-section becomes 
weighted away from small halos (despite arising preferen- 
tially in small galaxies) towards more massive halos, where 
forming structures and outflow-induced motions can yield 
higher velocities. This results in substantial amounts of 
super-gravitational motion that directly translates into wide 
DLA systems. 

HSR98 claimed agreement with observed DLA kinemat- 
ics by noting that their individual DLA simulation produced 
a velocity width of 60% of the virial velocity v v i r , and then 
integrating over the halo mass function assuming that the 
DLA cross-section scales as v^ iT . In our comparable no- wind 
simulation, this latter scaling is appropriate for larger halos, 
but it becomes much shallower for smaller ones, yielding 
much more DLA absorption at small cross sections. This re- 
sults in poor agreement with the observed kinematic s in the 
no-wind case. In another study, iMaller et al. I (|200lh found 
that gaseous disks of standard sizes were unable to reproduce 
DLA kinematics, and that more extended gaseous distribu- 
tions were required. While they mostly focused on testing 
extended disks models, they noted that outflows could po- 
tentially have a similar effect. Our momentum-driven wind 
case appears to produce just the required gas distribution. 
We note that disk sizes are generally too small without 
strong feedback and much hig her numerical resolut ion than 
we have here, but the results of lMaller et al. I (|200lf ) suggest 
that even optimally conserving angular momentum is not 
enough, and some supergravitational motion is required. Fi- 
nally, while outflows produce greater velocity widths, they 
also produce less cross-section overall. We hypothesize that 
such models will therefore produce more Lyman Lim i t Sys - 
tems (LLS; 10 17 < N m < 1 20 - 3 cm~ 2 ): iDave et all l|201Ch 
showed that winds do produce more strong Lya forest sys- 
tems (Nm > I0 14 cm~ 2 ) at low-z, but did not have suf- 
ficient statistics to explore LLS. CDM models have tra- 
ditionally also had difficulty producing enough LLS (e.g. 
iHernquist et al.l 1 1996h . though the strength of the discrep- 
ancy has been difficult to quantify owing to uncertainties in 
measuring column densities in this regime. 

More broadly, we note that momentum-driven winds 
have enjoyed substantial success in reproducing a wide range 
of observations of early galaxies and intergalactic gas, par- 
ticularly related to cosmic chemical evolution. On the other 



hand, this investigation of DLAs focuses solely on the neu- 
tral gas; metals are merely used as a tracer of dense Hi. 
Combining these various results leads us to conclude that 
outflows must move not only metals but also substantial 
amounts of mass around the cosmos. It may be possible to 
construct a model that enriches the IGM by preferentially 
expelling highly enriched gas from galaxies which could still 
reproduce the metal-line kinematics, but the strong low- 
ionization metal absorption generally indicates a substan- 
tial neutral hydrogen component, if for nothing else than 
to shield it from being sent into higher ionization states. 
Furthermore, the high mass outflow rates (typically com- 
parable to or greater than the mass forming into stars; 
lOppenheimer fc Dave I [2008 ?) suggest that a large amount 
of ISM gas must be entrained in these outflows, and hence 
the outflow metallicity cannot be substantially greater than 
the ambient ISM. Hence it is difficult to envision a sce- 
nario where the metals are being ejected highly preferen- 
tially compared to the mass. The fact that a single outflow 
model where vast amounts of gas are expelled at the typical 
ISM metallicity reproduces all these observations within a 
cosmological context is quite remarkable. Nevertheless, the 
lack of possibly important physics in our current simulations 
precludes any definitive constraints on outflows from DLA 
kinematics. 

DLA kinematics have long been an oddity that did not 
seem to fit neatly into our current hierarchical view of galaxy 
formation. Our study suggests that the answer is not in 
modifications to the hierarchical view, but rather to modi- 
fications in the detailed processes of how galaxies form. In 
a sense, this is a more exciting possibility, as it allows us 
to move towards constraining such processes from a com- 
pletely new and different perspective. The work presented 
here paints a hopeful picture that DLA kinematics may fit 
into a broader understanding of galactic outflows. However, 
our current form of modeling outflows is relatively crude, 
and the present analysis did not address the wider range of 
DLA observations such as metallicities and redshift evolu- 
tion that could provide interesting constraints. Furthermore, 
we did not include detailed radiative transfer to more ac- 
curately predict the distribution of neutral gas. Hence this 
work must be considered as a preliminary step that merely 
highlights the importance of a new physical process in yield- 
ing DLA properties as observed; we plan to examine a more 
comprehensive suite of DLA properties in the future. This 
will help us obtain a more complete picture for how DLAs 
relate to the evolution of galaxies in a hierarchical universe. 
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